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■ Abstract. The accurate simulation of collapsed objects requires that 
^4 \ a huge range of spatial scales be well resolved if the result is not to be 

■ contaminated by numerically induced fragmentation. In this context, 'in- 
', sufficient resolution' means comparable to the local instability wavelength 

■ in the system. I define a minimum resolution criterion for simulations of 
, disk systems, using the critical wavelength appropriate for disks. 

t-H ■ 

> 1 The Criteria in 2d and 3d 

o : 

, A condition on the minimum resolution to ensure the collapse in 3d is of 

^^•O ■ physical rather than numerical origin was defined by Truelove et al. 1997 using 

the ratio of the local grid resolution, 5x, and local Jeans wavelength, Aj: 

^: J = Yj where Aj = U;J • (1) 

9^: 

Here c s is the sound speed, p is the volume density and G is the gravitational 
+3 , constant. An analogous condition can be obtained for rotationally supported 

(i.e. disk) systems using the local Toomre wavelength, Aj-, at which the disk is 
neutrally stable to growing axisymmetric instabilities (see e.g. Lin & Lau 1979): 
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Here again, 5x is the grid size, while S is the local surface density. 

To obtain a form more useful in particle based numerical methods (e.g. 
SPH), Bate & Burkert (1997) using the Jeans mass to define an analogous cri- 
terion for the maximum resolvable volume density. We can obtain a similar 
quantity in 2d by defining a 'Toomre mass' and a maximum resolvable surface 
density. The two criteria are: 
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and S max = ( t« ——. $ • (3) 



N Teso is the number of particles (of mass m p ) required to resolve the Jeans mass. 

In 3d, Truelove et al. 1997 found emperically that J < 1/4 was sufficient to 
suppress numerical collapse. Bate & Burkert (1997) arrived at a similar result 
for particle simulations: numerical instability could be suppressed by resolving 
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Figure 1. Azimuth averaged (dotted line), maximum (solid) and crit- 
ical (red solid) surface densities for the Nelson et al. (1998) simulation 
modeled with 8000 (left) and 397000 (right) particles 



the local Jeans mass with N Teso > 2A^ nei gi 1 (~ 100) neighbor particles. I will 
define the same critical values for 2d as in 3d, i.e. T < 1/4 and ~ 100 particles 
per Toomre mass, saving a more exact specification for later work. 

Examples from the Literature 

Nelson et al. (1998) modeled the evolution of self gravitating disks with 
masses between 0.05 and 1.0 times the mass of the central star, using SPH. 
From the suite of models in that work, I choose as an example the model with 
disk dimensions Md/M*=0.2 and Q m \ n =1.5, resolved with 7997 particles, im- 
plementing an isothermal equation of state. Using physical dimensions, the disk 
has a radius i? D =50 AU, mass M D = 0.2M* and M* = O.5M . 

Figure 1 shows the surface densities produced at resolutions of 8000 (the 
original) and 397000 particles. At lower resolution, filamentary spiral arms de- 
velop that later collapse into multiple clumps, as in the Nelson et al. work, but 
such surface densities clearly violate criterion of eq. 3. The higher resolution 
model did not produce any clumps, even after evolving for three times as long 
(about 1500 yr). We originally concluded that although our physical model was 
insufficent to correctly model clump formation, if it was to occur at all it would 
be most likely between 10 and 40 AU. The statement that the physical model 
in Nelson et al. was insufficient is certainly valid (see e.g. Nelson et al. 2000). 
However, and clump formation was definitely invalid: it was due to failure of 
the Toomre criterion and not to any physical process, however modeled. 

Armitage &: Hanson (1999) showed collapse occured in disks in which a lMj 
point mass was embedded in a self-gravitating 0.1M© disk with an r _3//2 density 
profile and no initial gap, resolved with 60000 SPH particles. The minimum 
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Figure 2. Azimuth averaged (dotted line), maximum (solid) and crit- 
ical (red solid) surface densities for the Armitage & Hanson (1999) 
simulation. 



stability of the disk using the Toomre parameter, was Qmin — 1.5, as in Nelson 
et al. . Temperature was defined using a locally isothermal equation of state. 

Figure 2 shows a reproduction of this model. As in the original, the planet 
both accretes a large amount of disk matter and generates large amplitude spiral 
structures that began to collapse after ~ 20 — 25yr. Also as in the original, 
the initial radial positions correspond closely to the outer and inner Lindblad 
resonances for the lowest order spiral patterns, since the patterns reach their 
largest amplitudes there. These collapsed objects are clearly numerical artifacts 
and the author's conclusion that if a planet could grow sufficiently rapidly (to 
4-5Mj) sufficiently early enough in the disk's evolution, it would trigger further 
planet formation in the disk, is not supported. 

In the most recent of a series of papers studying collapse in disks, Boss 
(2002) models a disk in 3d (reflected through the midplane), between 4 and 
20 AU, with 100 x 23 x 256 grid zones asymmetrically distributed primarily near 
the disk midplane. He includes an ideal gas equation of state and a radiative 
cooling model. He evolves a disk for 359 yr, when the simulation begins to 
violate the Jeans criterion due to the onset of clump formation. 

Figure 3 shows the critical wavelength at 359 yr defined using the Jeans 
criterion (Boss: personal communication) and using eq. 2. The values from the 
two criteria coincide at the location of the clump. Elsewhere, the Toomre wave- 
length exceeds the Jeans wavelength. Thus the evolution of the Boss simulation 
leading up to its termination passes this test, and the conclusions made from it 
are valid numerically. 

Concluding Questions and Remaining Work 
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Figure 3. The failure of the Jeans criterion (solid line) was used by 
Boss (2002) to define the end of his simulation. The Toomre crite- 
rion applicable for circumstellar disks (heavy solid line) shows a near 
identical critical wavelength at this resolution. 



What is the proportionality factor for the 2d criterion? So far, 
I have only assumed that the spatial resolution required in the 3d case would 
carry over in 2d as well, but this may not be correct. 

Which critical wavelength criterion is more applicable: Toomre 
or Jeans? The universe is (spatially) a 3d place, so a criterion based on only 
two coordinates is necessarily approximate. On the other hand, the geometry of 
a disk is highly flattened. The assumption made in the Jeans analysis-of an infi- 
nite (in all directions) homogenous medium-is not satisfied and may incorrectly 
capture the physical behavior of the system. Further, correspondence between 
the two criteria in the forming clump of Boss's simulation may be coincidental. 
The two criteria have differing functional dependences on the temperature. 

Except for the simulation discussed in figure 3, generously shared prior to its 
publication by Alan Boss, the computations reported here were performed using 
the UK Astrophysical Fluids Facility (UKAFF). I also gratefully acknowledge 
financial support from UKAFF. 
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